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Using the Ginzburg-Landau theory for two-band superconductors, we determine the surface energy 
cTs between coexisting normal and superconducting states at the thermodynamic critical magnetic 
field. Close to the transition temperature, where the Ginzburg-Landau theory is applicable, we 
demonstrate that the two-band problem maps onto an effective single band problem. While the 
order parameters of the two bands may have different amplitudes in the homogeneous bulk, near Tc 
the Josephson-like coupling between the bands leads to the same spatial dependence of both order 
parameters near the interface. This finding puts into question the possibility of intermediate, so 
called type-1.5 superconductivity, in the regime where the Ginzburg-Landau theory applies. 
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I. INTRODUCTION 

Depending on the behavior in external magnetic fields, 
superconductors are classified as type-I or type-IL In 
type I superconductors, the surface energy density Og be- 
tween regions of finite and zero order parameters, coex- 
isting at the thermodynamic critical field He, is positiveii 
In type-II superconductors this energy is negative and a 
homogeneous superconducting state is no longer stable, 
leading to the formation of a vortex lattice.— Within the 
Ginzburg-Landau (GL) theorjii for one-band supercon- 
ductors the interface energy per unit area, 



T(K) 



4^('^)' 



(1) 



is determined by the value of the thermodynamic critical 
field, -ffcj the magnetic penetration depth. A, and the 
dimensionless function, T (k) , that depends on the GL 
parameter k = A/^, the ratio of the penetration depth 
and the superconducting coherence length. Properties of 
T (k) are discussed, e.g. in ReL 0. In the regimes of 
extreme type-I and type-II superconductivity 



T(^) = 



(V2- 



if K < 1 
if K > 1 



(2) 



We have evaluated this function numerically and the 
result is shown in Fig. 1. The transition between type-I 
and type-II behavior occurs for k = 2^^/^, where T (k) 
changes sign. 

Fermi surfaces in many superconductors may con- 
sist of two or more well separated sheets with differ- 
ent energy gaps4i^ Evidence for two energy gaps was 
obtained in high-purity superconducting Nb, Ta, and 
V,^ and Nb-doped SrTiOa.^ Recently, tunnelingSi^iS and 
point contact^^'^^ spectroscopies, as well as heat capacity 
measurements^^^^'5 for MgBa^""— give clear evidence for 
two-band superconductivity with gaps Ai ~ 0.7meV and 
A2 ~ 2.5meV (for recent reviews, see Refs. [2l| andl22|). 
Other systems that have been discussed as two-band sys- 
tems are RNiaBzC with R=Lu,Y^i24 2H-NbSe2^ and 
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Figure 1: (Color online) The function, T («:), of the one-band 
problem calculated numerically (full points) compared to the- 
oretical limit of Eq. [2] (dashed lines). The inset shows an 
enlargement of the large k, domain. 



the recently discovered FeAs superconductors.^^ In all 
cases the amplitude of the superconducting gap is differ- 
ent for different sheets of the Fermi surface. 

Motivated by the study of these multi-band supercon- 
ductors, the term type-1.5 superconductivity has been 
coined^l to emphasize the possibility of a state that is 
intermediate between the two regimes. Specifically, one 
considers two-component or two-band systems with or- 
der parameters ^1 (r) and ^2 (r) that have qualitatively 
different spatial dependence, with different respective co- 
herence lengths ^1 and ^2- The existence of these two 
length scales emerges from the assumption that one can 
neglect the Josephson type coupling between two order 
parameters. The regime where one expects novel behav- 
ior is obviously the limit ^1 ^ A ^ ^2- Then one order 
parameter component might behave as a type-I super- 
conductor while the other follows the type-II behavior. 
Consequences of such behavior were discussed in Ref . |28| 
where it was concluded that properties emerge that fall 
outside the usual type-I/type-II dichotomy. For exam- 



pie, the emergence of "vortex molecules" and of an inho- 
mogeneous state comprising a mixture of domains of a 
two-component Meissner state and vortex clusters were 
proposed. In Ref. '29 the surface energy for such a sys- 
tem was analyzed with the conclusion that T (k) must 
be replaced by a function that depends on several di- 
mensionless quantities, in particular on the ratio ^1/^2- 
Changing ^1/^2 at fixed penetration depth was shown to 
yield a sign change of dg . 

Obviously, even in multiband superconductors, the 
sign of the surface energy is either positive or nega- 
tive. Thus, it seems more appropriate to discuss the 
physics that was investigated in Ref. [23 within the GL 
approach,— as interesting modifications of type-II super- 
conductivity. More importantly, it is crucial to analyze 
what exactly happens in a multiband superconductor in 
the vicinity of the transition temperature, with 



r-(T, -T)/Te«l, 



(3) 



in the regime where the GL approach is valid (ignoring, 
as usual, critical fluctuations). 

One of the key features of the two-band GL model is 
the Josephson like coupling between the two bands. 



/, (r) = -r; (*J (r) *2 (r) + % (r) *i (r)) , 



(4) 



in the expansion of the GL free energy density. Refs. [2^ 
and|29| analyze the limit ?? = 0, but assume that both or- 
der parameters, while uncoupled, order at the same tem- 
perature. The more realistic regime is clearly the one 
where the common transition temperature is the conse- 
quence of a finite order parameter coupling 77. 

In this paper we determine the surface energy Us of 
a two-band GL model including the coupling, Eq.Q, 
between the bands, i.e. we consider rj ^ 0. We find 
that in the regime r ^ 1, where the GL theory provides 
the correct mean field description, Eq.((T]) continues to be 
the correct expression for the interface energy with same 
function T (k), which implies that the surface energy con- 
tinues to change sign at k = 2~^/^. The multi-component 
nature of the order parameter enters the GL k through 
the values of A and 



e = (Cr- 



<^r'^'- 



(5) 



A detailed definition of A and ^i in the two-band prob- 
lem is presented below. We also find that, while the 
order parameters may have different amplitudes in the 
homogeneous bulk, ^i_o and ^2,0, close to the transition 
temperature, i.e. for small r, they have the same spatial 
dependence near the interface. In particular: 



^2(2) * 



2,0 



0(r) 



(6) 



i.e. the coupling enforces the same spatial dependence 
for both components. ^1 (r) and 'i'2 (r) vary in space 
on the single length scale ^ of Eq. ([S]) . Close to a supercon- 
ducting transition it is then sufficient to introduce only 



one order parameter to characterize the symmetry bro- 
ken state. An exception is the case where two completely 
uncoupled order parameters are accidentally degenerate, 
i.e. both components accidentally have the exact same 
transition temperatures Tc, while they have, at the same 
time, different coherence lengths. This is the scenario 
considered in Refs. [2^ and [23. We stress, that our re- 
sults do not preclude novel type-II behavior that may 
occur deeper in the ordered state away of the GL do- 
main. This is however beyond the limit of applicability 
of the GL theory. In the next section we present our 
analysis. We summarize our findings in section 3. 



II. TWO-BAND GINZBURG-LANDAU THEORY 



We start with the free energy, 

F = J f{r)d^ 



(7) 



of a two-band superconductor. .F is a functional of the 
pairing wave functions VPi and \l/2 of the two components 
or bands and of the vector potential A associated with 
the magnetic field B = V x A. The free energy density, 
/ (r), relative to the zero field normal state value, is: 



/(r) = /i(r) + /2(r) + /,(r) + 






(8) 



Here the fj (r), with j = 1, 2, are the GL expansions of 
the two bands: 



fj 



Jjl*j 



^'« 



I*, 
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2m* 
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(9) 
and fc (r) is the coupling term given in Eq. (jj]) . Here 
bj > and the bands' effective masses are tuj. The physi- 
cal values of the order parameter and vector potential are 
determined via SF/6'^i = SF/SA^ = 0. In principle ad- 
ditional coupling terms such as (*i'I'2)^ , V\E'* • V*2, 
etc. are allowed. For clean multiband systems, a weak 
coupling expansion yields that the coefficients of such 
terms vanish due to momentum conservation^^. In ad- 
dition, even if present, such terms are sub-leading close 
to the transition temperature point when compared to 
/eWofEq.O. 

We first discuss the homogeneous, zero field solution. 
Ignoring the inter-band coupling, fc, one finds, as usual, 
^i.o iv = 0) = \f^a~Jh^ for a, < and f.^o (77 = 0) = 
for fli > 0. In the general case of /c 7^ 0, however, the 
common critical temperature T^ is not equal to either of 
Tc^j and there is no reason that both coefficients a^ (T) 
change sign at the same temperature. At T^ the smallest 
eigenvalue of the matrix of the homogeneous quadratic 
terms in / (r) vanishes. In our problem, this eigenvalue 
is 



r- = i (ai+a2- V(ai-a2)' + V); (10) 



it vanishes for rf^ — ai (Tc) 02 (Tc). Thus, it must hold 
that ai^2 [Tc) > 0, as r_ would be negative if one of the 
two ai is smaller or equal to zero. Thus, close to T^ 
tti > and the interband coupling enhances the transi- 
tion temperature compared to the largest of the Tcj for 
the 77 = limit. 

To proceed, we introduce the dimensionless ratio 



t = 



1]- 



aia-z Tc-T 
ex 



aia2 



Tr 



(11) 



that vanishes at Tc (see also the Appendix). It is conve- 
nient to eliminate 77^ = (1 -I- t) 0102 in favor of t. Thus, 
small t naturally corresponds to finite interband coupling 
77. For small t we have r_ ~ —^0102/ (ai -I- 02) and the 
smallest eigenvalue changes sign at t = 0. 

The free energy minimization of the homogeneous 
problem for 77 7^ leads to the system 




Figure 2: (Color online) Schematic representation of the in- 
terface between the normal and the superconducting state. 



ai*i -t-61*? -77*2 = 
02*2 + ^2*2 -?7*i = 



(12) Using Eq.([T4]), we obtain 



which is readily reduced to a fourth order equation for ^^ 
that can be solved using known formulas for the roots of 
a quartic equation. One can simplify the problem by rec- 
ognizing that the GL theory is only valid in the vicinity of 
the transition temperature, t <C 1. The homogeneous or- 
der parameters can easily be determined to leading order 
in t: 



^1.0 



*2,0 



uit with ui 

U2t with U2 



a|ai 



0261 -I- afb2 

a\a2 
a^bi + a\b2 ' 



(13) 
(14) 



where the subscript is to denote the zero-field solution. 
Hence, the temperature dependence of the order parame- 
ters is as expected: 



*,o cxicx 



Tr~T 



Tr 



(15) 



We stress that within GL theory there is no reason to 
go to terms of higher orders in t. Of course, away from 
Tc corrections can be significant, in particular for small 
77, but those effects require a microscopic approach based 
on Bogoliubov-de Gennes equations^i. 

Close to Tc, we can also determine the thermodynamic 
critical field by imposing / (He) = 0: 



A^^ = 



47re* 0102 a2/m\ + ai/m^ 
(? a^bi + a\b2 



(18) 



It is now straightforward to set up the formalism to de- 
termine the interface energy. 



III. THE INTERFACE ENERGY 

In evaluating the surface energy we follow closely the 
classical approach that was used for the single band 
problem. "'^ Consider the interface between superconduct- 
ing and normal half-spaces at the plane 2 = 0. The field 
H is applied along the x axis parrallel to the interface and 
equal to He to ensure coexistence of two phases. Then the 
magnetic induction has only one component B.^ ~ B (z) 
and the vector potential can be chosen as Ay — —A (z), 
as shown in Fig. 2, yielding 



B{z)^A'{z). 



(19) 



We can choose real order parameters that vary along 
the z-direction. The Gibbs free energy per unit area that 
is minimum in a given applied field reads: 



4tt 



dz B (z) . 



(20) 
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(16) 



One can formally define the one-band penetration depth 

v-2 _ ^_„*2,T,2 / r^*„2\ 



as Aj ^ = 47re*^^^g/ {m*c^). Since the additive super- 
fluid density is proportional to A~^, the actual London 
penetration depth is 



K^+^2^- 



(17) 



Far away from the interface we have on one side the nor- 
mal state with \l/i (z -^ —00) -^ and B {z —> — cx)) ^• 
He, while on the other side for z — > 00 we have the 
homogeneous superconducting state with S — > and 

It is convenient to introduce dimensionless quantities: 



V' 



\I/2 
J 



B 



V2Hc' 



and s 



z 
A' 



(21) 



which imply that the dimensionless vector potential a 
A/ (\/2-ffcA) . The surface energy is then given by 



and effective parameter k given by 



H 

CTs = X^^[tpi, ^2, a], 



(22) 



where S is a functional that must be minimized with 
respect to the ipi, ^p2 and a to yield as. After simple 
algebra we obtain: 

S = y f l^ (Vi>2, a) + ^ ^,7 Vf + (a' - 2-1/2) 'J ^^ 

(23) 



where 



^(■01, ■02, a) 



V^i + V'l - 2VT+lVl02 , u 



r. 



l__u 4 ^ nltPl + KJi )? 



2 „2 



a'. (24) 



We use here the following notations: ijj'^ — dipi/ds, b = 
a' — da/ds, u = wi, as given in Eq. p^ . and 



K, = — , With ti — . 



(25) 



Minimization of S gives a system of coupled differential 
equations for tpi, a: 



nf' 2 90, 



IdV 
2~da' 



(26) 



(27) 



Multiplying Eq.lgni) by -0^, Eq.d!?]) by a' and summing, 
the first integral of this system is obtained: 

Yl '^i^^? +a'^ -V {^1 ,^2,a)= const. (28) 



The peculiar term in our analysis is the first one in 
V {ipi ,ip2,o,) of Eq. ((24)) that seems to be singular as t —> 
0. Expanding for small t, we have: 



^? + 0i - 2ym0iV2 (v-i - ^2)' 



i 



i 



V'iV'2. (29) 



Thus, close to the transition temperature we must have 
■01 = '4'2. Introducing ^ (s) = '0i (s) = "02(5), which 
is equivalent to Eq.®, one obtains the surface energy 
functional in the form: 

Y.= f ds (Vo (^, a) + K- V^ + [a - 2-^/2) j , (30) 

with 



Voiij,a) = -i^^ +^4^V'a' 



,-2 _ ,.-2 



(32) 



This is an exact form of the functional for the standard 
one-band surface enrgy problemjii^ 

It is worth noting that Ki enter the surface energy only 
through the combination k of Eq.(|32|l. In particular, this 
leads to Eq. ^ for the correlation length of the two band 
problem with k ~ X/^. Thus, the interface problem is 
identical to the one of a single band system, leading to 
Eq.lJT]) with the same function T (k). 

These conclusions are supported by numerical min- 
imization of S [^1 ,^2,o]- We discretized the interval 
s — [0, 2L] to N equidistant steps {sj = 2jL/N) and 
minimized S with respect to "01 (si), ^2 (si) and a{si) 
subject to boundary conditions a (0) = 2-^/^, a {2L) — 0, 
and i/'i (0) = and tpi (2L) = tpi^o- The homogeneous 
bulk solutions V'i.o approach the value 4'i,o = 1 for i — >■ 0. 
Finally, since in the limit 2L -^ 00 the interface position 
is arbitrary, a.t z = L we assumed ipi (L) = i, which 
centers the interface position in the large n limit. 

Our results for N — 400 are shown in Figs. 3-6. In 
comparing Fig. 3 and Fig. 4, as well as Fig. 5 and Fig. 
6, we show that the order parameters do indeed approach 
the behavior with identical spatial variation, as given by 
Eq-®, as the critical temperature is approached. 

In Figs. 3 and 4 we focus on the most nontrivial 
limit with ki = 0.45 < 2"^/^ < ^2 = 5. Naively, 
one could expect "01 to change on distances of the order 
^1 > -\/2A (type-I behavior), while ^2 < V2X suggests 
type-II behavior of "02 • Contrary to this expectation we 
find that both order parameters are strongly coupled by 
the Josephson energy and have increasingly similar spa- 
tial variation as i — )■ 0. As we will see below, the interface 
energy for this set of parameters is positive and the sys- 
tem behaves as a type-I superconductor as k in Eq. p2p 
is dominated by the smallest of the two Ki . 

In Figs. 5 and 6 we show the behavior for ki — 3 and 
K2 — 4, corresponding indeed to a type-II superconduc- 
tor with negative interface energy (see below for explicit 
values). Again, both order parameters follow the same 
spatial dependence and behave according to Eq.® as t 
decreases. 

In addition, as t decreases, the value of the mini- 
mized functional of Eq. (P5|) approaches the value of the 
function T (k) of the single band problem with k de- 
termined by Ea. p2p . This can explicitly be seen in 
the numerical results of "Smin and T (k) correspond- 
ing to Figs. 3-6. The effective one-band solution with 
K = 0.448, thus corresponding to Fig. 3 and Fig. 4, 
gives T (k) — 0.479, which differs from the numerical re- 
sult of Fig. 3, E„„ = 0.530, by ~ 11%. This difference 
decreases to ~ 2% for a smaller value of t as shown in 
Fig. 4, for which Smm — 0.488. The numerical solu- 
tions shown in Fig. 5 obtained for t = 0.2 correspond to 
Smin — —0.275 whereas the effective one-band k = 2.4 
(31) yields T (k) = -0.47: hence E™i„ and T differ by - 42%. 
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Figure 3: (Color online) (a) The reduced field and order para- 
meters obtained numerically by minimizing the interface en- 
ergy functional, E, for k\ = 0.45, K2 = 5, u = 0.6, t = 0.07. 
The reduced order parameter ■i/'i is shown by a solid line, the 
dashed line is i/;2. (b) The close-up of the order parameters 
for 6.5 < z/\ < 7.5. 



Figure 4: (Color online) (a) The reduced field and order para- 
meters obtained numerically by minimizing the interface en- 
ergy functional, E, for aci = 0.45, K2 = 5, u = 0.6, t = 0.01. 
The reduced order parameter tpi is shown by a solid line, the 
dashed line is 'tp2- (b) The close-up of the order parameters 
for 6.5 < z/X < 7.5. 



Again, by decreasing the value of t to 0.01, we find re- 
sults shown in Fig. 6 corresponding to Smm = —0.459, 
now only by '-^ 2% different from T. 



IV. SUMMARY 

In summary, for a two-band superconductor we ana- 
lyzed the energy of the interface between regions of a 
finite order parameter and zero order parameters, co- 
existing at the thermodynamic critical field H^- If one 
includes the interband Josephson coupling between the 
bands, i.e. the leading allowed interaction between the 
Cooper-pair wave function \E'i and \I'2, both order para- 
meters vary close to the transition temperature on iden- 
tical length scales. Thus, despite the fact that both order 



parameters may have very different amplitudes, they vary 

1/9 

on the same characteristic length scale (^;j~^ + ^2'^) , 
where the ^i are the typical length scales where the gra- 
dient (or kinetic) energies in the GL functional become 
comparable to the bulk condensation energy. We stress 
that ^1,2 are just auxiliary quantities and only ^ is a mea- 
surable physical length. An important implication of this 
result is that the surface energy is determined by a single 
GL parameter k = X/^ in a way identical to the single 
band case. 

Thus, there is no room left for so-called type- 1.5 su- 
perconductivity in the GL regime close to Tc- Of course, 
our analysis cannot rule out the possibility of interest- 
ing novel physics due to distinct charactetristic length 
scales S.i deep in the superconducting state. Such possi- 
bility then requires an approach within the microscopic 
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Figure 5: (Color online) (a) The reduced field and order pa- 
rameters obtained numerically by minimizing the interface 
energy functional, E, for ki = 3, K2 = 4, u = 0.6, t — 0.2. 
The reduced order parameter xpi is shown by a solid line, the 
dashed line is i/j2. (b) The close-up of the order parameters 
for 6.5 < z/\ < 7.5. 



Figure 6: (Color online) (a) The reduced field and order pa- 
rameters obtained numerically by minimizing the interface 
energy functional, E, for ki = 3, K2 = 4, u — 0.6, t — 0.01. 
The reduced order parameter ■i/'i is shown by a solid line, the 
dashed line is ip^. (b) The close-up of the order parameters 
for 6.5 < z/\ < 7.5. 



framework of the Bogoliubov-de Gennes or Gor'kov for- 
malisms. In the recent review by Brandt and Das,"^^ sit- 
uations are described which do not fit to a rigid type-I- 
type-II dichotomy. Close to Tc, however, only one rel- 
evant superconducting order parameter exists and the 
phenomenology of the transition between type-I and 
type-II superconductivity is unchanged by the multiband 
character of the system. 
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Within a weak coupling BCS theory it is possible to 
derive the parameters of the GL expansion, Eqs. ([2l)-([9|), 
in terms of the microscopic densities of states and pairing 



interactions)^ 



with effective coupling constant 



6,; = 



V 



87r2T2 "*' 

Nf , 



In 



TrTr 



detA 



12- 



(Al) 



Here Nf is the densities of states at the Fermi level 
per one spin, Ui = NF,i/NF are relative densities of 
states on two bands, A^- — NFVij are interaction con- 
stants proportional to the symmetric matrix Vij respon- 
sible for the Cooper pairing, ^^^ r = {Tc — T) jT^, and 
(A~^) = A22/detA etc. are elements of the matrix in- 
verse to \ij . 

The transition temperature follows from the condition 
\o-\o-'2\r^Q = ^ that leads to 



Tc = W£) exp ( -1/Aj 



(A2) 



A = 2niri2detA [riiAii + n2A22 

- \J (niXii - n2\22f + 4nin2Af2 



(A3) 



It is now straightforward to express the variable t in 
Eg. dlip in terms of t close to the transition temperature, 
which shows 



det A-^ (niAii - n2A22) + 4?iin2Af. 

'^12 
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